AD-A1 62  029 
UNCLASSIFIED 


RELIA8ILITV  ASSESSMENT  FOR  SYSTEMS  SUBJECT  TO  121 

MAINTENANCE  AND  REPAIRIU)  CITV  COLL  NEU  YORK  DEPT  OF 
MATHEMATICS  M  BROUN  AUG  85  AFOSR-TR-85-1884 
AFOSR-84-0095  F2G  1424  NL 


PROCUREMENT  instrument  identification  NUMIER 


A  NAME  OF  FUNOINQ/SPONSORING 
ORGANIZATION 

AFOSR 


k-  ADDRESS  (City.  State  an*  ZIP  Co* I 

Bolling  AFB,  D.C.  20332-6448 


AFOSR-84-0095 


10.  SOURCE  OF  FUNDING  NOS. 


PROJECT 

NO. 


WORK  UNIT 
NO. 


13.  PERSONAL  AUTHOR  (SI 

Mark  Brown 


13a  tvfe  of  refort 

Annual 


13A  Time  COVERED  1 1«.  DATE  OF  RSFORT  (Yt.,  Mo..  Da yt  16.  FAQE  COUNT 

fromI  Apr  84  tc80  Jun  8E  August  1985  7 


COSAT  I  CODES 


GROUP 


SUS.  OR. 


IS  SUBJECT  TERMS  /Conttnm  oa  raaan*  if  MMBary  ana  Mr nttfy  by  Mock  aannkari 

^Reliability  modeling,  age-weighted  distributions, 
repairable  systems. 


IS.  ASST R ACT  rCoallMM  on  raaaiw  1/ memory  an*  Mmfl/y  ky  Mock  namkarj 

*>The  main  thrust  of  the  research  under  this  grant  is  to  develop  the  methodology  to  deal 
with  more  realistic  models  for  the  quantitative  reliability  assessment  and  analysis  of 
complex  systems.  Inequalities  for  increasing  failure  rate  distributions  have  been  derived. 
A  new  theory  of  age-weighted  distributions  is  being  developed.  Work  is  beginning  on 
exploiting  properties  of  a  pratial  ordering  of  distributions.  DTIC 

tf&ELECTEf| 

V  DEC  O'!  1985  Js  i 


m  fiiE 


30i  OISTRISUTION/AVAILASILITV  OF  AMTRACT 


21.  AMTRACT  SECURITY  CLASSIFICATION 


UNCLAMlF  IE  D/UN  LIMITED  Q  SAME  AS  RFT.  E  OTIC  USERS  □ 


Unclassified 


JflFOSR-TR 


,  O  ;  . 


Interim  Scientific  Report  to  the  Air  Force 
Office  of  Scientific  Research 


Grant  AFOSR-84-0095 


Reliability  Assessment  for  Systems  Subject 
to  Maintenance  and  Repair 


Department  of  Mathematics 
The  City  College,  CUNY 
New  York,  NY  10031 


Mark  Brown,  Principal  Investigator 


August,  1985 


y  ■  .  v  • 


-  •  V/.V-1  •• 


mm 


m 


5‘W 


AIR  tORCS  0?P1  C" 

*071  Cp.  . . 

7hi'  >  ;  ,1  ■ 


Grant  AFOSR-84-0095  was  awarded  for  the  period  4/1/84  -  3/31/85. 

It  was  subsequently  extended  to  6/30/85  without  further  funding.  A 
renewal  proposal  covering  7/1/85-6/30/86  was  submitted  and  approved. 
This  report  will  discuss  the  research  accomplishments  so  far  achieved 
during  the  grant  period. 

The  main  thrust  of  the  research  performed  under  the  grant  is  the 
development  of  methodology  to  deal  with  more  realistic  models  for  the 
quantitative  reliability  assessment  and  analysis  of  complex  systems. 

This  methodology  is  potentially  applicable  to  the  programs  of  the  U.S. 
Air  Force  in  particular  and  the  Department  of  Defense  in  general.  In 
addition  the  results  obtained  are  of  value  in  a  variety  of  other  applied 
areas,  and  in  probability  and  statistical  theory. 


1)  Inequalities  for  IFR  Distributions. 


A  technical  report  Inequalities  for  distributions  with  increasing 
failure  rate"  was  issued  in  December  1984.  The  paper  will  appear  in  a 
special  volume  in  celebration  of  the  65th  birthday  of  Herbert  Solomon. 

Several  new  inequalities  for  IFR  (distributions)  are  derived.  Included 
is  a  two  sided  bound  on  the  renewal  function  as  well  as  the  exact  sup  norm 
distance  between  the  renewal  function  and  its  asymptotic  linear  approximation 
This  is  one  of  several  papers  the  author  has  written  on  the  subject 
of  exponential  approximation.  The  importance  of  this  study  is  that  quite 
often  the  first  two  moments  of  a  distribution  .are  known  along  with  the 
reliability  class.  The  above  approximations  give  uniform  two  sided  bounds 
for  the  cdf  which  are  excellent  when  the  parameter  o  is  small.  As  the 


distributions  of  interest  are  often  intractable,  the  ability  to  do  such 
approximation  is  invaluable.  In  particular  it  allows  one  to  closely 
approximate  the  distribution  of  the  time  to  first  failure  for  repairable 
systems. 


2)  Age-weighted  distributions. 

Under  the  current  grant  I  have  been  developing  a  new  theory  which 
looks  interesting  and  potentially  useful.  I  plan  to  write  a  technical 
report  in  the  near  future  summarizing  the  current  state  of  development. 

Given  a  distribuiton  F  on  [0,°°)  with  finite  mean  p  define 
AF,  the  age-weighted  distribution  of  F  by: 

d(AF) (t)  »  ~  dF(t)  . 

Age-weighted  distributions  arise  in  many  different  areas.  Two  promi¬ 
nent  examples  are: 

(i)  If  F  is  the  interarrival  time  distribution  for  a  renewal 
process,  then  AF  is  the  distribution  of  the  length  of  the  interval 
covering  the  point  0  in  a  stationary  renewal  process  on  the  whole 
real  line  with  interarrival  distribution  F.  The  fact  that  F  and  AF 
have  different  distributions,  with  AF  being  larger,  is  the  content  of 
the  waiting  time  paradox  discussed  by  Feller  (1966). 

(ii)  Suppose  that  'v-  F^  and  X£  ^  with  and 

mutually  absolutely  continuous.  Let  f^  and  f0  be  the  densities 
of  F^  and  ?2  with  respect  to  \  (a  measure  dominating  both  F^ 
and  ^2^’  Then  the  Neyman-Pearson  test  of  vs  F£  is  based  on 


2  2 

the  ratio  Define  Yi  “  f-  and  Y?  =  “  0^)  >  the  Neyman- 

Pearson  statistic  when  F^  and  F2,  respectively,  are  true.  Then 
Y2  =  AY^,  i.e.  the  distribution  of  the  Neyman-Pearson  likelihood 
ratio  under  F2  is  the  age-weighted  distribution  of  its  distribution 
under  F^. 


Below  are  some  of  the  results  I  have  obtained  during  the  current 
grant  period  for  age-weighted  distributions. 

Let  X  'V  F  and  T  ^  AF.  Then: 


(1) 


ET  =  (EX) 


-1 


More  generally,  let  X^,...,Xq  be  i.i.d.  as  X  and  independent 


of  T.  Define  S  =  -  X . .  Then: 

n  i»l  i 


Accesion  For 


(2) 


E(S  +T)"-1  =  [  (n+DEX]"1 
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The  relationship  between  X  and  T  ^  is  symmetric,  i.e. 


(3) 
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.-1 


-1 
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T=  (AX)  and  X  =  (AT  ) 


Defining  c2  =  1-(EXEX_1)  1  and  c2 
from  (l)-(3)  that: 
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1-(ETET  )  it  follows 


(A) 


c2  *VarT/ET2 


(5) 


cT  =  Var  X/EX  . 


Finally  as  a  consequence  of  (l)-(5),  let  g  be  a  convex  function 

3 

with  w  =  sup(x  g"(x))  <  ®.  Then: 
g 

(6)  0  <  Eg(X)-g(EX)  <j  (EX‘1-(EX)"1)wg  . 

The  import  of  (6)  is  that  X  may  have  a  distribution  with  large 
or  infinite  variance.  As  a  result  the  standard  delta  method  for  approxi¬ 
mating  Eg(X)  may  not  be  applicable  or  useful.  However,  if  g  is 
rapidly  decreasing,  the  mean  of  g(X)  should  not  be  extremely  sensitive 
to  the  tail  of  the  distribution  of  X.  Relationship  (6)  gives  an  alter¬ 
native  to  the  delta  method  in  which  the  variance  (which  measures  the 

2  2  -1 
departure  of  EX  from  (EX)  ) ,  is  replaced  by  the  departure  of  EX 

from  (EX)  This  appears  to  be  a  fundamental  idea,  well  worth  deeper 

study. 

3)  Partial  ordering  of  distributions. 

In  previous  papers  I  obtained  some  interesting  results  for  IMRL  and 
DFR  distributions  by  exploiting  a  partial  ordering  between  F  (IMRL  or 
DFR)  and  its  stationary  renewal  distribution  G. 

I  am  currently  abstracting  this  partial  ordering  and  developing  its 
properties.  It  looks  like  a  promising  new  tool.  Below  are  some  of  my 
results. 


Let  Fi>F2  probability  distributions  on  [0,®).  Define 
_>  if  F1(t)/F2(t)  is  increasing  in  t  for  F2(t)  >  0.  It 
turns  out  that  this  partial  ordering  is  equivalent  to  each  of: 

(i)  Xx  =  min(X2,Z)  where  X1  v  F^,  X2  “v  F2  and  Z  is  a 
(possibly  improper)  random  variable,  independent  of  X2- 

(ii)  _>  G2  is  the  sense  of  monotone  likelihood  ratio,  where 
G^  is  the  stationary  renewal  distribution  corresponding  to  F^, 
i  -  1,2. 

(iii)  Let  X^  ^  F^,  X2  ^  F2«  Then  for  all  t  2.  0,  the  conditional 
distribution  of  X^-t  given  X^  >  t  is  stochastically  larger  than  the 
conditional  distribution  of  X2~t  given  X2  >  t. 

The  main  results  I  have  derived  for  this  partial  ordering  are: 

(a)  Let  F^  >  F2,  with  F^  continuous  and  H^Ct)  *  -ln(F^(t)). 

Then : 

*» 

suplF.  (B)-F-(B)  [  <_  1  -  F,(t)dH.  (t) 

BeS  J0 

where  B  is  the  collection  of  Borel  subsets  of  [0,®). 

(b)  Let  X^,  X2  be  independent  with  X^  v  F^,  X2  v  F2  and 
F^  <_  F2»  Let  h(x,y)  be  a  function  of  two  variables  with 
h(x,y)-h(y ,x)  increasing  in  x  ^  y.  Then: 


(7) 


Eh(X. ,X9)  >  Eh(X-,,X- )  . 
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